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We investigate the viability of coupled waveguides as basic units of quantum circuits. In particular, 
we study the dynamics of entanglement for the single photon state, and single mode squeezed vacuum 
state. We further consider the case of entangled inputs in terms of the two mode mode squeezed 
vacuum states and the two photon NOON state. We present explicit analytical results for the 
measure of entanglement in terms of the logarithmic negativity. We also address the effect of loss on 
entanglement dynamics of waveguide modes. Our results indicate that the waveguide structures are 
reasonably robust against the effect of loss and thus quite appropriate for quantum architectures as 
well as for the study of coherent phenomena like random walks. Our analysis is based on realistic 
structures used currently. 
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I. INTRODUCTION 

Discrete optical systems like coupled waveguides are 
known to be extremely efficient in manipulating the flow 
of light and have been investigated extensively in the last 
two decades |2|, y, |j] . Many key quantum effects like 
quantum interference, entanglement and quantum walk 
has been investigated in these systems d, [f| 0, HJl] • For 
example, using coherent beam Peretes et al. [7[ have 
observed quantum walk effects in a system consisting 
of large number of waveguides. In another experiment, 
Bromberg et al. investigated the quantum correlations in 
GaAs waveguide arrays [8( using two-photon input states. 
In particular, they considered both the separable and en- 
tangled two photon state and observed various features 
associated with quantum interference. In addition, the 
coupled waveguides arrays have been used to study the 
discrete analogue of the Talbot effect [l(|. Note that, 
entanglement between the waveguide modes and behav- 
ior of nonclassical light in coupled wave guid es have also 
attracted a great deal of interest [ill fl2lTT3j . In a recent 
experiment Politi et. al. [l4| have shown how a CNOT 
gate can be implemented on a single Silicon chip using 
coupled silica waveguides, thus showing possible appli- 
cation of waveguides in quantum computation. They 
also observed two photon interference in these coupled 
waveguides. In a following experiment [Tol [l6| coupled 
silica waveguides was used to generate a multimode in- 
terferometer on an integrated chip. It was further shown 
that these interferometers can be used to generate arbi- 
trary quantum circuits. They also showed that two and 
four photon entangled states similar to NOON states [l7| 
can be generated on the silicon chip. All these studies 
have hence given a new impetus in the field of quantum 
information processing and quantum optics with wave- 
guides. Note that the entanglement between waveguide 
modes is at the heart of many of these experiments. In 
particular, for effective use of these waveguide circuits 
in quantum computation and communication tasks sus- 
tainability of generated entanglement is very important 
. [l9j. In light of this, it is imperative to study entan- 



glement in waveguides using quantitative measures for 
entanglement. This is the main purpose of the present 
study. Moreover, in practice the waveguides are not com- 
pletely lossless. Thus an immediate question of interest 
would be how does this loss affects the entanglement in 
the waveguide modes ? It is well known that entangle- 
ment is quite susceptible to decoherence [2(| and thus 
the above question bears immense interest in context to 
quantum information processing using waveguides. Fur- 
ther it is important to understand the role of loss in co- 
herent phenomena like quantum random walk 0, [2l| . 

In this paper we investigate these in a simple sys- 
tem of two single mode waveguides, which are coupled 
through the overlap of evanescent fields. This simple sys- 
tem serves as a unit or the basic element for construct- 
ing a quantum circuit (22J. The input light to the cou- 
pled waveguide system is usually produced by a paramet- 
ric down-conversion process at high and low gain which 
produces important nonclassical states of light like the 
squeezed and the single photon states respectively. Thus 
the input is quite naturally a squeezed state specially at 
high gain. Behavior of photon number states such as the 
single photon state and the NOON state have also been 
investigated in these systems [3, [H, H3 ■ We thus con- 
sider a variety of nonclassical input states like squeezed 
states and photon number states which have been exten- 
sively investigated in couple waveguide system and study 
their respective entanglement dynamics. We quantify the 
evolution of entanglement in terms of logarithmic neg- 
ativity and present explicit analytical results for both 
squeezed and number state inputs. We further investi- 
gate the question of possible effects of loss on the en- 
tanglement dynamics in waveguides by considering lossy 
waveguide modes. We find that in this case, for both 
number state inputs as well as squeezed state inputs, en- 
tanglement shows considerable robustness against loss. 

The organization of the paper is as follows: In Sec. II, 
we describe the model and derive analytical result for the 
field modes of the coupled waveguide system. In Sec. Ill, 
we study the evolution of entanglement for two classes of 
photon number states, (A) separable single photon state 
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1 1,1) and (B) entangled two-photon NOON state. We 
quantify the degree of entanglement of these states by 
using the logarithmic negativity. In sec. IV , we then 
study the time evolution of entanglement by evaluating 
the logarithmic negativity for two classes of squeezed in- 
put states (A) separable two mode squeezed state and 
(B) entangled two mode squeezed state. The effect of 
loss in waveguides on the entanglement dynamics is then 
discussed in Sec. V. Finally we summarize our results in 
section VI with a future outlook. 



description for the nonclassical light can be studied by 
quantizing the field amplitudes as has been done in the 
work of Lai et. al. [251 ]. Following an approach simi- 
lar to that developed by Lai et. al., we can write the 
corresponding quantum mechanical Hamiltonian for the 
coupled waveguide as 

H = hw(a)a + tfb) + hJ{a ] b + tfa) , (1) 



II. THE MODEL 

We consider a system with two single mode wave- 
guides, coupled through nearest-neighbor interaction as 
shown in Fig. 1. Let a and b be the field operators for the 
modes in each waveguide. These obey bosonic commuta- 
tion relations [a,a^] = 1; (a — > b). The Hamiltonian de- 
scribing the evanescent coupling between the waveguide 
mode in such a system of two coupled wa veg uides can 
be derived using the coupled mode theory [23, [24| . The 
coupling among the waveguides is incorporated in this 
framework by treating it as a perturbation to the mode 
amplitudes. It is assumed that the presence of the second 
waveguide perturbs the medium outside the first wave- 
guide. This creates a source of polarization outside the 
first waveguide, which thereby leads to modification of 
the amplitude of the mode in it. Further, the amplitude 
of the modes in each waveguide is assumed to be a slowly 
varying function of the propagation distance. Moreover, 
in this perturbative approach the coupling does not ef- 
fect the propagation constant or transverse spatial dis- 
tribution of the waveguide modes. The field of the first 
waveguide has a similar effect on the second waveguide. 
Under these assumptions, the field mode of the composite 
structure are governed by the Helmholtz equation which 
gives two coupled first order differential equations which 
can be solved to obtain the time evolution of field modes 
in the coupled waveguide structure. The corresponding 



where the first two terms correspond to the free energy of 
the waveguide modes and the last two terms account for 
the evanescent coupling between the waveguide modes 
with J as the coupling strength. The coupling J depend 
on the distance between the waveguides. The input to 
the coupled waveguide system can be in a separable or 
an entangled state. Let 7 be the loss rates of the modes 
a and b. The loss 7 arises from the loss in the material 
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FIG. 1: (Color online) Schematic diagram of a coupled wave- 
guide system. The parameter J gives the coupling between 
the waveguide modes and 7 is the loss rate. 

of the waveguide. Table I below gives the experimental 
values of coupling parameter J and loss 7 for different 
waveguide systems. 



Waveguide Type 


Coupling parameter J (sec x ) 


Loss 7 (sec : ) 


7 /J 


Lithium Niobate (LiNb0 3 ) 


1.83 x 10 10 - 4.92 x 10 10 


3 x 10 9 


1/7-1/20 


AlGaAs 


2.46 x 10 11 


2.7 x 10 10 


1/10 


Silica 


1.53 x 10 11 


3 x 10 9 


1/50 



TABLE I: Approximate values of some of the parameters used in waveguide structures [26l. I27I. |28|] . The loss, usually quoted 
in dB/cm, for different waveguides is converted to frequency units used in this paper by using the formula, 10 Log( F ° ut ) = 

10 Log(e~ 2l ^ c ), where Pi„ is the input power, P ou t is the power after traveling unit length. 

I 



As known the silica waveguides have very little intrinsic Nevertheless the loss is to be included as this could 
loss and should be preferable in many applications. be detrimental in long propagation for example in 
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the study of quantum random walks. Since the two 
waveguides are identical, we have taken the loss rate 
of both the modes to be the same. We can model the 
loss in waveguides in the framework of system-reservoir 
interaction well known in quantum optics and is given by, 

Cp = —^{a^ap — 2apa) J rpa)a) 

~(b r bp-2bptf + ptfb) , (2) 

where p is the density operator corresponding to the 
system consisting of fields in the modes a and b. The 
dynamical evolution of any measurable (O) in the cou- 
pled waveguide system is then governed by the quantum- 
Louiville equation of motion given by, 

P=-~[H,p] + Cp (3) 

where (O) = Tr{Op}, the commutator gives the uni- 
tary time evolution of the system under the influence 
of coupling and the last term account for the loss. Note 
that in absence of loss (lossless waveguides) the time evo- 
lution of the field operators can be evaluated using the 
Heisenberg equation of motion and is given by, 

a(t) = a(0) cos( Jt) - ib{0) sin( Jt) 

b(t) = 6(0) cos( Jt) -m(0)sin(Jt). (4) 

Next we will study the entanglement characteristics of 
photon number and squeezed input states as they propa- 
gate through the waveguides. To keep the analysis simple 
in the next few sections we consider the case of lossless 
waveguide modes (7 = 0). We defer the discussion of loss 
on entanglement to Sec V. 



III. EVOLUTION OF ENTANGLEMENT FOR 
NON-GAUSSIAN INPUT STATES 

In this section we study the dynamics of entanglement 
for photon number input state. We quantify the entan- 
glement of the system by studying the time evolution for 
the logarithmic negativity (26|, [30, [M|. For a bipartite 
system described by the density matrix p the logarith- 
mic negativity is 

Bv(t)=log a ||p T ||, 

|| p T ||= (2N(p) + 1) , (5) 

where p T is the partial transpose of p and the symbol || 
denotes the trace norm. Also N(p) is the absolute value 
of the sum of all the negative eigenvalues of the par- 
tial transpose of p. The log negativity is a non-negative 
quantity and a non-zero value of Ej^ would mean that 
the state is entangled. 



A. Separable photon number state as an input 

We first consider the case when there is no loss and hence 
we set 7 = 0. We assume that the input is in a separable 
state. Further, for studying the entanglement dynamics 
for photon number states we first consider the case of a 
single photon input in each waveguide. Thus the initial 
state is 

|#))> = |1,1>. (6) 

Using Eq. (|4|) we can show that a single photon input 
state given by \ip(0)} evolves into a state : 

m)) -ai|2,0) + A|l, 1) + *i|0,2) . (7) 

The coefficients ax, Pi and S\ are given by : 

ai = -isin(2J<)/V2, 
/3 X = cos(2Ji), 

5i = -isin(2Jt)/V2 . (8) 

The density matrix corresponding to the state in (JT]) can 
be written as : 

p = \m)(m\ 

= M 2 |2,0)(2,0| + lAplMXMl + N 2 |0,2)<0,2| 
+a 1 %\2, 0)<1, 1| + *i/3f|0, 2)(1, 1| + ai Sl\2, 0)(0, 2| 
+Pia$\l, 1)(2,0| 1)(0,2| + 5 ia t\Q, 2)(2,0| . 

(9) 

Using Eq. ([5]) and the above equation, we can show that 
the log negativity Ej^ is given by: 

E M =log 2 (l + 2JV(p)), 

= log 2 (l + 2|(a 1 /3 1 + aitfi + 5i/3i)|) . (10) 




0=Jt/7T-> 

FIG. 2: Time evolution of log negatively for a single photon 
input state 

In Fig. ([2]) we show the time evolution of E^r for the single 
photon input state |1,1). We would like to emphasize 
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that the values of 9 studied here are very similar to the 
ones employed in the recent experiments [H, [l4[ . At time 
t = 0, we begin with a separable input state and thus 
the value of log negativity is Erf = 0. The entanglement 
quantified by the log negativity increases with time and 
attains a maximum value of 1.58 for 9 ~ 0.15. In this 
case the single photon state evolves into a maximally 
entangled state given by: \ip m ) = e~ l7T / 2 (\2, 0) + |0, 2)) + 
|1, 1)/V3. Further, for 9 = 1/4, we get an analog of 
the well known Hong-Ou-Mandel interference 32]. Note 
that in this case the logarithmic negativity Erf attains a 
value of 1 which is less than the corresponding value of 
Erf for the maximally entangled state \ip m ). In addition, 
for 9 = 1/2, we find that Erf vanishes and the state at 
this point is e^\l, 1). At later times, we see a periodic 
behavior which can be attributed to the inter-waveguide 
coupling J. We next consider the case where we have two 
photons in one waveguide and none in the other input. 
Thus the initial state can be written as : 

|p(Q)> = |2,Q>. (11) 

Again using Eq. ((4| we find that the |v(0)) evolves into 
a state given by : 

\<p(t)) ->a 2 |2,0)+/?2|l,l> + <fa|0,2) . (12) 
The coefficients ot2, 02 and 62 are given by : 
a2 = cos(Ji) 2 , 
P2 = ~V2icos(Jt) sin(Ji), 
82 = -sin(Jf) 2 . (13) 

Using a similar procedure as discussed above we can eval- 
uate the log negativity Erf for the state in Eq. (|12[). We 
show the result for the log negativity in Fig. ([3|). In this 
case we find that the log negativity increases and attains 
a maximum value of 1.54. After reaching the maximum 
value the log negativity decreases and eventually becomes 
equal to zero. Thus the state becomes disentangled at 
this point of time. At later times we see a periodic be- 
havior and the system gets entangled and disentangled 
periodically. Clearly the entanglement dynamics of the 
states (JSJ) and (fTTj) are different. Unlike the earlier case 
for the 1 1,1) input state, we don't see any interference 
effects in this case 0. 

B. Entangled photon number state as an input 

Next we consider the entangled state prepared in a two 
photon NOON state [13] as our initial state : 

i, ( o), . m±m . (14) 

As shown in the black curve of Fig. (TJ|, the value of 
Erf at time t = is equal to 1 which indicates entangle- 
ment. The log negativity Erf in Fig. shows a behavior 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

0=Jt/7T-> 

FIG. 3: The behavior of log negatively for the state (|lip as 
function of 9 = Jt/n. 




0=Jt/7T-» 

FIG. 4: The behavior of log negatively as function of 8 = 
Jt/n. The black curve shows the result for two photon NOON 
state while the red curve shows the result(_Ejv — 1) for the four 
photon NOON state 



that is similar to the result for the |1, 1) state shown in 
Fig. ((2]). Also note the shift of 7r/4 between the results 
in Figs. (f5]) and ([?]). As in the case of single photon in- 
put state |1, 1), the initial state evolves into a maximally 
entangled state corresponding to a value of Erf which is 
equal to 1.58. In addition, for 9 = 1/2, we again see a sig- 
nature of quantum interference such that the probability 
of gettingthe single photons at each of the output port 
vanishes [8[ . The logarithmic negativity Erf at this point 
is equal to 1. At later times the entanglement shows an 
oscillatory behavior and the system gets periodically en- 
tangled and disentangled. For an arbitrary NOON state 
: \ip in ) = (\N,Q) + \0,N))/V2 the state at time t is 
\ipout) = (52Pk\k,N — k)) where 0u can be written in 
terms of the binomial coefficient (see Eq. (|A3jl ). For com- 
pleteness, we give the details for the calculation of time 
evolution of \ip% n ) in Appendix A. The density matrix cor- 
responding to the state \ip ou t) can be written as : p ou t = 
^2 PkPm\ki N — k) (m, N — m\. Taking the partial trans- 
pose of p ou t, we get p T out = Y / (3k(3* l \k,N-m) (m,N-k\. 
Further it can be proved that (p^ ut ) is a diagonal matrix 

and the eigenvalues of (p£ ut ) is of the form: |/?fc| 2 |/3 m | 2 . 
Thus the negative eigenvalues of p^ ut are of the form 
|/3fc||/3 m | (k 7^ to) and the log negativity Erf can be writ- 
ten as: 
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where a, (3 and \i are 2x2 matrices given by, 



E N =tog 2 (l + 2JV(p)), 

= log 2 (l + 2^|^||/3 m |) . 

We can use the above equation to study entanglement 
dynamics for the N photon NOON state. The red curve 
in Fig. 4 shows the result for the four photon NOON 
state. As earlier, the value of £jv at time t = is equal 
to 1 which indicates entanglement. The curve for four 
photon NOON state also shows quantum interference ef- 
fect. Further, the logarithmic negativity never becomes 
zero in this case and hence the initially entangled state 
remains entangled for later times. 

IV. EVOLUTION OF ENTANGLEMENT FOR 
GAUSSIAN INPUT STATES 

A. Separable two mode squeezed state as an input 

We next study the generation and evolution of entan- 
glement for the case of squeezed input states. For this 
purpose we first consider a separable squeezed input state 
coupled to the modes a and b of the waveguide given by, 

10 = IC«> ® 10,}; (15) 

where |Ca)(lO>)) are single mode squeezed states defined 
as, 



ICa) 



exp(-{a 



12 



2 })|0>; (o-&). 



(16) 



where r is taken to be real. It is well known that a two 
mode squeezed state like |£) can be completely character- 
ized by its first and second statistical moments given by 
the first moment : ((xi), (pi), (£2), {P2}) and the covari- 
ance matrix a. The squeezed vacuum state falls under 
the class of Gaussian states. It is to be noted that evolu- 
tion of Gaussian states has been studied for many differ- 
ent model Hamiltonians [H, [H, HH, [36| . We focus on the 
practical case of propagation of light produced by a down 
converter in coupled waveguides which currently are used 
in quantum architectures and quantum random walks. 
Note that since the first statistical moments can be arbi- 
trarily adjusted by local unitary operations, it does not 
affect any property related to entanglement or mixedness 
and thus the behavior of the covariance matrix a is all 
important for the study of entanglement. The measure of 
entanglement for a Gaussian state is best characterized 
by the logarithmic negativity _Eyv, a quantity evaluated 
in terms of the symplectic eigenvalues of the covariance 
matrix a (37l . |38| . The elements of the covariance matrix 
<7 are given in terms of conjugate observables, x and p in 
the form, 



a fi 



(17) 



\ 2 



V = 



1 x 1 pi+p 1 x 1 \ 
* 2 ' 



> X 2 P2+P2X 2 \ 
i 2 / 



/ XXX2+X2X! ' 
\ 2 / 



/ X2P1+P1X2 \ 

\ 2 / 



(pi) 



I X2P2+P2X2 \ 
\ 2 / 



(Pi) 



(18) 



(19) 



I X1P2+P2X1 > 

\ 2 1 



I P1P2+P2P1 ' 
\ 2 - 



(20) 



Here X\,Xi and pi,P2 are given in terms of the nor- 
malized bosonic annihilation (creation) operators a(a^), 
b{w) associated with the modes a and b respectively, 



Pi 



(a + a f ) 
V2 ' 

(o - at) 
V2t ' 



X2 



P2 



(fr + fot) 
y/2i 



(21) 



The observables, Xj,pj satisfy the cannonical commuta- 
tion relation [a^,^] = iSkj- The condition for entangle- 
ment of a Gaussian state like \Q is derived from the PPT 
criterion [38[ , according to which the smallest symplectic 
eigenvalue v < of the transpose of matrix a should satisfy, 



where v < is defined as, 



v < = min[i/_|_, V— 
and v± is given by, 



v± 



\ 



A(cr) ± J A(cr) 2 - 4Detcr 



(22) 



(23) 



(24) 



where A(cr) = A(ct) = Det(a) + Det(/3) - 2Det(/x). Thus 
according to the condition (I22p when i>< > 1/2 a Gaus- 
sian state become separable. The corresponding quantifi- 
cation of entanglement is given by the logarithmic nega- 
tivity Em [H S3 defined as, 



E N {t) = max[0,-ln{2i><(f)}]; 



(25) 



which constitute an upper bound to the distillable en- 
tanglement of any Gaussian state [39l |. On evaluating 
the covariance matrix a for the state (|15p for 7 = (no 
loss), using equation ([3]), (JH) and (|2"Tj) we find, 



(26) 



c 




e 




; p = 
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e 
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FIG. 5: Plot of the time dependent logarithmic negativity 
E^j- for the state \C). Here amount of squeezing is taken to be 
r = 0.9. 



where d,e,c are given by 



c = -{cosh(2r) + sinh(2r)cos(2Jt)}; 
d = i{cosh(2r) - sinh(2r) cos(2J<)}; 



e = --sinh(2r)sin(2Jt). 



(27) 



The corresponding symplectic eigenvalues v± are then 
given by 



(28) 



v± = Vcd ± e 



One can clearly see from equations (|25[) . (|27|) and (|28|) 
the dependence of logarithmic negativity -Eyv' on coupling 
strength J between the waveguides and the squeezing pa- 
rameter r. In figure 5. we plot the logarithmic negativ- 
ity as a function of scaled time, 8 = Jt for the state \(). 
Here t is related to the length I of the waveguide and its 
refractive index n by t — nl/v, v being the velocity of 
light. We see from figure (5) that as |C) is separable at 
t = 0, Ej\f = initially but as Jt increases, it oscillates 
periodically between a non-zero and zero value. Thus 
the initially separable state \Q becomes periodically en- 
tangled and disentangled as its propagates through the 
waveguide. We attribute this periodic generation of en- 
tanglement to the coupling J among the waveguides. We 
further find that v < = 1/2 at certain points along the 

waveguide given by 28 = (k+l)n, k = 0, 1, 2, 3, Note 

that at this points Ej^ vanishes and |£) becomes separa- 
ble. At all other points the state \Q ^ \Q a ) ® |Cb). We see 
that E^f is maximum and has a value equal to the amount 
of squeezing 2r at the points given by 28 = (k + 1)tt/2. 
Hence at this points the initial seperable state |£) be- 
comes maximally entangled and is given by, 



|C) =exp{e i,r r(a t 6 t + a6)}|00) 



(29) 



B. Entangled two mode squeezed state as an input 

Let us now study the dynamical evolution of a two mode 
squeezed state |£) as an input to the waveguide, 



10 =exp[r(a t 6 t -a&)]|00) 



(30) 



As before we consider r to be real. To quantify the entan- 
glement of the state |£) we need to evaluate the logarith- 
mic negativity E^f. Thus we first evaluate the covariance 
matrix a for the state |£) using equations ([3]) with 7 = 0, 
flU) and (dU). We find a to be 



/ 9 h 

g f o -h 

h / g 
-h g f 



(31) 



where /, g and h are given by, 



1 



/ = -cosh(2r) 
g = -- sinh(2r) sin(2 Jt) 
h = i sinh(2r) cos(2 Jt); 



(32) 



0.5 



2.5 



FIG. 6: Time evolution of logarithmic negativity EV for the 
initial entangled state Here the squeezing is taken to be 
r = 0.9. We see that is 7r/4 out of phase to that for the 
state \(). 



The corresponding symplectic eigenvalues v± is then 
given by, 



= y/(f + g)(f-g)±h; 



(33) 



The logarithmic negativity E^f can then be evaluated us- 
ing equations (f2"3")l . (|23|) and (|3"3")) . From equations (|3"2"|) 
and (|33|) the dependence of £V on the squeezing r and 
the coupling J between the waveguides is clearly visi- 
ble . From equation (|33|) we find that Ej^r — i.e. 
entanglement become zero when, 28 = (k + l)n/2 as 
then v < = 1/2 and thus the initially entangled state |£) 



becomes separable, i.e |£) = exp{|e l7r (a t2 



2 )}|0) ® 



exp{§e i7r (^ 2 + 6 2 )}|0). 
evolution of E^/ for r 



I 'i figure (6) we plot the time 
0.9. We see that entanglement 
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oscillates periodically between zero and non-zero values. 
We further find that in this case the oscillations in E^ 
is 7r/4 out of phase to that for the initial separable state 
|C) . This oscillatory behavior of entanglement is as dis- 
cussed before, due to the coupling J among the wave- 
guides. Each time the states get separable the presence 
of coupling leads to interaction among the modes of the 
waveguides and creats back the entanglement. We see 
from the figure that logarithmic negativity Etf reaches 
maximum at later times at the points 29 = (k + 1)tt and 
is equal to 2r. Thus at this points the state |£) regains 
its initial form given by equation (|30p . 



V. LOSSY WAVEGUIDES 

In this section we study the entanglement dynamics 
of lossy waveguides (7 7^ 0). The loss 7 arises from 
the loss in the material of the waveguide. In this case 
the dynamical evolution of the waveguide modes is gov- 
erned by the full quantum-Louivillc equation ([3]). We 
next consider the cases of both photon number state and 
squeezed states at the input of the waveguide and discuss 
the influence of the loss on their respective entanglement 
evolution. 



A. Effect of Loss on non-Gaussian Entanglement 

As discussed above, we first study the effect of loss 
on the entanglement dynamics of the waveguide modes 



for photon number input states. For this purpose we 
consider a single photon input state |1,1) as the ini- 
tial state. In this case we can analytically solve the 
quantum-Louiville equation described in (J3]). To pro- 
ceed further, we work in the interaction picture such 
that the density matrix in the interaction picture is 
p(t) = e iJt ( atb+bU )p(t)e- iJt( - aH+bU \ In the interaction 
picture we can write Eq. ([3|) as 



dm 

dt 



7 , , , 

- — (a'ap — 2apa' + pa' a) 

-^(¥b p - 2bptf + pftb) , (34) 



where a and b are given by 



a(t) = acos(Jt) — ibsin(Jt) 
b(t) = bcos(Jt) — iasin(Jt). 



(35) 



Using the above equation, we can rewrite Eq. (|34|) as 



em = _7 
2 

_7 
2 



dt 



(a'ap — 2apa) + pa'a) 
{tfbp-2bptf + ptfb) . (36) 



J 



For the separable input state |1, 1), the solution for the density matrix (|36p can be written as 4l| : 



p(t) = e - 47t {(e 27 *-l) 2 |0,0)(0,0| + (e 2 ^-l)|l,0>(l,0| + (e 2 ^-l)|0,l)(0,l| + |l,l)(l,l|}. (37) 

I 



Further, we can write p(t) in terms p(t) using the follow- 
ing equation: 

p(t) = e -^t(^b+b*a)~^yJt(aH+^a) ( 38 ) 

The above equation gives the time evolution of the den- 
sity matrix corresponding to the single photon state 
|1, 1). Following a similar approach as discussed in Sec. 
Ill, we can evaluate the log negativity for the lossy wave- 
guide case. But the resulting expressions are lengthy 
and do not exhibit a simple structure. Thus we only 
give the numerical results for the lossy waveguide case. 
In Fig. ([7]), we show the decay of entanglement, as a 
function of scaled time for the state ©■ Note that the 
range of 7/ J values studied here are similar to the nu- 
merical values used in the experiments [26l . [27j ■ For ex- 
ample, the coupling parameter J for the lithium niobate 
waveguide lie between 1.83 x 10 10 sec -1 and 4.92 x 10 10 



sec -1 . The loss parameter for these waveguides is close 
to 3 x 10 9 sec -1 [20| which corresponds to a value of 7/J 
between 1/7 and 1/20. For AlGaAs waveguides the loss 
7 is close to 2.7 x 10 10 sec -1 [27| . The coupling param- 
eter J for these waveguides is about 2.46 x 10 11 sec" 1 . 
Thus the 7/J value for these waveguides is of the order of 
1/10. It is worth mentioning that the 7/J value for silica 
waveguides is significantly lower than the corresponding 
values for the lithium niobate and AlGaAs waveguides. 
This means that even a small loss would add up to a 
significant decoherence in these complex quantum sys- 
tems. From Fig. ([7]) we find that for the lossy waveguide 
case the entanglement between the waveguide modes de- 
crease with time. In addition, we find that increasing 
the value of 7/J makes the waveguide modes more frag- 
ile, as is evident from Fig. ([7]). However, we find that 
the decrease in entanglement is not substantial. Our re- 
sults indicate that the waveguide system can sustain the 



entanglement even for the higher decay rates. Thus the 
coupled waveguide system can be used as an efficient tool 
for the study of basic quantum optical effects. In addi- 
tion, the persistence of entanglement suggests that the 
coupled waveguide system can be used effectively for var- 
ious applications in quantum information processing [lH ] . 
For example, the single photon entanglement described 
here is a key step for the successful implementation of 
the CNOT gate [1J| . We also studied the behavior of log 
negativity for the entangled initial state (|14[) . In this case 
also we found that the entanglement quantified by Ejy 
shows a considerable robustness against the decoherence 
effect. 
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FIG. 7: (Color online) Time evolution of the logarithmic neg- 
ativity Em in presence of loss of the waveguide modes for the 
initial separable input state ((6)1 . The decay rates of the modes 
are given by 7/J = 0.1 (solid black), 7/J = 0.2 (broken black) 
and 7/J = 0.3 (red). 



B. Effect of Loss on Gaussian Entanglement 

For the input squeezed state \Q of equation (fl~5|) we 
find that elements of the covariance matrix er in presence 
of loss become dependent on the decay rate 7 and is given 
by, 



d e' 

d! e' 

e' J 

e' d! 



(39) 



where c , d! , e are given by 

1 
2 
1 

2 

1 



r = -{l + e~ 27t sinh 2 (r) + e - 27 *sinh(2r)cos(2Jt)}; 
{1 + e~ 27t sinh 2 (r) - e~ 27t sinh(2r) cos(2Jt)}; 
27 *sinh(2r)sin(2Jt). (40) 



d' 



The corresponding symplectic eigenvalue v of the co- 
variance matrix is then found to be, 



LU 
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FIG. 8: (Color online) Time evolution of the logarithmic neg- 
ativity Ej\r in presence of loss of the waveguide modes for 
the input state \C). The decay rates of the modes are given 
by 7/J ~ 0.1 (solid black), 7/J = 0.2 (broken black) and 
7/J = 0.3 (red). Here the squeezing is taken to be r — 0.9. 
The loss leads to new behavior in the entanglement. 



v± 



(41) 



On substituting equation (I4lj) in equations ([23]) and us- 
ing (f2l)]) we get the logarithmic negativity-for lossy wave- 
guides. To study the dependence of entanglement on loss 
of the waveguide modes we plot the logarithmic negativ- 
ity Erf for different decay rates 7/J in figure (8). As for 
the case of single photon states we focus on the range of 
8 important from the experiment point of view. We see 
new features in the entanglement dynamics as an effect 
of the loss. We see from figure (8) that in presence of 
loss the maximum value of entanglement for the state |£) 
reduces in comparison to the case of lossless waveguides. 
However it is important to note that this decrease is not 
substantial. We further find that with increase in de- 
cay rate, the entanglement maximum shifts but does not 
show considerable reduction (the maximum changes by 
only 0.4 as the decay rate becomes three times). Thus 
we sec that entanglement is quite robust against decoher- 
ence in this coupled waveguide systems. The robustness 
of entanglement dynamics is an artifact of coherent cou- 
pling among and the waveguide modes. This findings 
hence suggest that coupled waveguide can be used as an 
effective quantum circuit for use in quantum information 
computations. Further we see another new feature in 
entanglement in figure (8). We find that there exist an 
interval of 8 during which the state |£) remains separable. 
Note that in absence of loss the state \Q becomes sepa- 
rable momentarily and entanglement starts to build up 
instantaneously once it becomes zero (see figure 5.) Thus 
this feature that entanglement remains zero for certain 
interval of time arises solely due to loss. 
In figure (9) we plot the long time behavior for entan- 
glement of the state |£) with very small decay rate of 
7/J = 0.1 and squeezing parameter r = 0.9. We see 
that entanglement decays slowly with increasing 8 as the 
magnitude of E^f diminish successively with every os- 
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e = jt/n 



FIG. 9: Long time behavior of the logarithmic negativity Em 
in presence of loss of the waveguide modes for the initial sep- 
arable input state \(). 




= Jt/Jt 

FIG. 10: Time evolution of the logarithmic negativity £V in 
presence of loss of the waveguide modes for the initial entan- 
gled input state Here 7 is the decay rate of the modes 
and squeezing is taken to be r = 0.9. 



where f',g', h! are given by, 

/' = i+e- 2 ^sinh 2 (r) 

g' = -ie" 27 *sinh(2r)sin(2Ji) 

h' = ^e~ 27t sinh(2r) cos(2 Jt); (43) 

In this case we find that the symplectic eigenvalues u± 
are dependent on the decay rate of the waveguide modes 
and is given by, 

v ± = ^m + m- ± h! (44) 

where m±{t) = 1 — e _27 *[l — {cosh(2r) ± 
sinh(2r) sin(2J£)}]. The corresponding measure of 
entanglement given by the logarithmic negativity EV 
can then be calculated by using equation |44|) . (f23j) and 
(|25|) . In figure (10) we plot the logarithmic negativity 
Ej\f for the state |£) as a function of in presence of loss. 
We find similar behavior in the entanglement dynamics 
as seen earlier for the separable state \Q. We find in 
figure (10) that entanglement of the state |£) decrease 
slowly with increasing 9 for non-zero 7/ J. Thus as for 
the separable states, in case of initial entangled input 
states entanglement is found to be quite robust in the 
face of loss. In addition to this we also see in figure (10) 
periods of disentanglement appearing successively as 9 
increases. 

The loss in waveguides that we discussed in this section 
arises due to material properties like change in refractive 
index and absorption. On the other hand there can be 
decay of the waveguide modes in the form of leakage to 
its surrounding also. It should be noted that leakage is 
inherently different from the evanescent coupling as the 
former can arise due to scattering and refraction due to 
refractive index difference at the waveguide boundaries. 
Thus the analysis of this section is also valid when the 
leakage is important as for example is the case when one 
couples channel waveguides to slab waveguides 0, [43| . 



dilations. In addition periods of disentanglement arises 
repeatedly in its oscillations. We find that the length of 
this periods increases with increasing 9. It is worth men- 
tioning here that this kind of behavior has been predicted 
earlier for two qubit entanglement (42l |. 

Next we study the effect of the decay of waveguide 
mode on the entanglement dynamics of the initial entan- 
gled squeezed state |£) given in equation ([317)) . We find 
in this case the covariance matrix to be, 



/' g' h' 

g' f -h' 

h' /' g> 

-h' g' f 



(42) 



VI. CONCLUSION 

To conclude, we investigated the time evolution 
of entanglement in a coupled waveguide system. We 
quantified the degree of entanglement between the 
waveguide modes in terms of logarithmic negativity. 
We have given explicit analytical results for logarithmic 
negativity in case of initially separable single photon 
states and for separable as well as entangled squeezed 
states. We have also addressed the question of decoher- 
ence in coupled waveguide systems by considering loss 
of waveguide modes. For the lossy waveguides we found 
that the entanglement shows considerable robustness 
even for substantial loss. Note that our results are 
based on experimental parameters and thus should 
be relevant for applications of waveguides in quantum 
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information sciences. Our results serve as guide for 
experiments dealing with entanglement in waveguide 
structures. For efficient use of these waveguides, one 
should choose the waveguide parameter like 8 so that one 
is away from values where the entanglement is minimum. 



APPENDIX A: TIME EVOLUTION OF THE 
INITIAL NOON STATE (\N,0) + \0,N)/V2) 



I 

In this appendix we give the details of our calculation for \ipout) when the input state is given by 



\lpin) = 



(|tf,0> + |0,JV» _ ((a(0)) w + (6(0)n|0,0» 



\/2 y/2N\ 
Using Eq. (j4|) we can show that the input state given by \ipi n ) evolves into a state 



(Al) 



IV'o 



(( a (t))" + (&(t)mo,o)) 

V2~/V! 



(A2) 



where a(t) and b(t) are given by Eq. J4|. Using Eq. (j4| in the above equation, we get 



IVw) = (%2l3 k \k,N-k)) , 

0k = a k + OtN-h , 

a k = (C(^,fc)) 1/2 (cos(Jt)) fc (-isin(Ji)) JV - fc , (A3) 
where C(N, k) is the Binomial coefficient given by: C(N, k) = Nl/(N — fc)!fc! . 

This work was supported by NSF grant no CCF-0829860. 
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